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We propose a quantum theory of swimming for swimmers that are small relative 
to the coherence length of the medium. The quantum swimming equation is derived 
from known results on quantum pumps. For a one-dimensional Fermi gas at zero 
temperature we find that swimming is topological: The distance covered in one 
swimming stroke is quantized in half integer multiples of the Fermi wave length. 
Moreover, one can swim without dissipation. 

The theory of classical swimming studies how a cyclic change in the shape of a swimmer 
immersed in a fluid leads to a change in its location. The theory is both elegant and practical 



and has been applied to the swimming and flying of organisms, robots 



and microbots 



A classical medium may be viewed as a limiting case of quantum medium when the 
quantum coherence length is small compared to the size of the swimmer and interference is 
negligible. Quantum mechanics takes over once the swimmer is small enough. Our primary 
goal here is to present a theory of swimming in quantum media. We concentrate on swimming 
in a one dimensional ideal Fermi gas at low temperatures but our methods are more general 
and can be applied to higher dimensions. One dimension is distinguished in that interference 
effects are especially strong (and also leads to tractable and soluble models). 

The swimmer, a q-swimmer, is an object with internal degrees of freedom which, we 
assume, are slow and classical while the degrees of freedom of the medium are fast and 
quantum. Swimming is accomplished by the internal degrees of freedom, the controls, un- 
dergoing periodic cycles. An example of a q-swimmer might be a molecule immersed in a 
quantum gas. The control of the internal configuration might be either external (through 
the application of external fields) or internal (due to autonomous dynamics). 

Natural q-swimmers and artifical q-bots may be regarded as quantum machines that 
propagate by pumping the quantum particles of the ambient medium. This point of view 



reflects an intimate connection between quantum pumps 
a role below. 



and q-swimming which will play 
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Consider a simple model of a swimmer made of n disconnected spheres of radii aj im- 
mersed in either a classical or quantum medium. The swimmer can control the n — 1 relative 
distances ii = Xi — X„, i = 1, . . . , n — 1 between the centers Xi of the spheres and the radii 
aj. Allow the swimmer to change adiabatically the control parameters aj and ii. The notion 
of adiabaticity means that the velocities, e.g. Xj, aj, are small compared with the charac- 
teristic velocities of the particles in the medium. By linear response, the force on the j-th 
sphere is given by 

fj = -Yl ^^^^^ ~ ^^^^^ + ^j^- 

k k k 

Here Fji = —F^ are internal forces acting between i and j spheres and rjjk, fjk are coefficients, 
which, a-priori, depend on the state of the swimmer (i.e. the relative distances li and the 
radii of the spheres aj) and the nature of the medium. 

To derive an explicit form of the swimming equation one first needs to make a choice how 
to designate the position of the swimmer, X. For a swimmer made of n disconnected pieces 
it is convenient to pick X = Xn, the coordinate of one of the components. The total force 
X]j fj acting on an adiabatic swimmer must vanish. (This follows from the fact that the 
friction forces are of first order in the adiabaticity while the acceleration is second order.) 
This constraint determines the swimming equation 

n— 1 n 

—T] dX = rji Mi + Vj daj , (2) 
i=i j=i 

where rj^ = J2]=i Vjk , V = X]j=i Vj ; ^k = Yl]=i ^jk- Swimming is manifestly geometric being 
independent of the (time) parametrization of the swimming stroke. The notation dX stresses 
that the position of the swimmer will, in general, not integrate to a function on the space 
of controls: X will not return to its original values when the controls undergo a cycle. 

Eq. (j2I) does not determine the coefficients rji, uj and rj. In this sense, the swimming 
equation, although general, is incomplete. We start by deriving a new formula for the 
quantum friction ?7 [7] . t] expresses the friction acting on an idle q-swimmer while it is being 
dragged at small velocity through the ambient quantum medium. The formula is in the 
spirit of Landauer formula ^,0]; It is expressed in terms of the scattering data. As usual 
in the Landauer setting we assume a one dimensional system. For the sake of simplicity, we 
focus on the two channel case. 

The most general on-shell scattering matrix of a one dimensional time-reversal invariant 
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FIG. 1: A scatterer with two scattering channels: {r,t) and {r',t') are the reflection and transmis- 
sion amplitudes 

scatterer anchored at the origin, can be parameterized as follows jl^ : 

i sin 9 e~*° cos 9 \ 
cos 9 i sin 9 I 

The reflection r, r' and transmission t, t' amplitudes are functions of three independent real 
parameters: a, 7,6'. 

The scattering matrix of a scatterer located at X is related to the scattering matrix 
located at the origin, Sq, by: 

S{X) = e'^^/^^oe-'''''/^ (4) 
where the on-shell momentum matrix P is deflned by 

(5) 




p{E) = \/2mE for electron gas and p{E) = E/c for photons. 

A dragged scatterer may be thought of as a pump. If the velocity of dragging is small 
compared with the characteristic velocity of the scattered particles, the theory of adiabatic 
pumps can be applied. In particular, the rate of momentum transfer to the ambient medium 
is |ll|: 

Vit) = --L f dEp'{E)TrEU{E,t)p), (6) 

- .... . e!:;; .... Q. . „ a .o.„ o..e„ 

scattering matrix and its time derivative. Here and after S* denotes the Hermitian conjugate 
of Tte denotes the trace on the scattering channels at fixed energy E and p{E) gives the 
occupation of the (scattering) states at energy E. If the ambient quantum gas is at thermal 
equilibrium, p{E) = (^e^^^~^^) ± is the Fermi-Dirac or Bose-Einstein distribution. 
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In dragging a scatterer, the time dependence of 5* comes solely from the change of position, 
X. From Eq. (jl)) S = X[P, S]S* . Eq. (jH)) determines the force / on the swimmer: f = V. 

Define, as usual, the friction coefficient, r/, by / = —rjX. Combined with Eq. (jnj we get 
a Landauer type formula for the quantum friction 

V = -^jdEp\E) TvE ( [P, S] [P, Sr). (7) 

At thermal equilibrium p{E) is a decreasing function of the energy. This implies that t] is 
non- negative. 

For a Fermi gas at zero temperature p'{E) = ~6{E — Ep). In the two channels case 
one then has rjp = -^p^^Ep) Ir^Ep)]"^. The friction depends only on the momentum and 
reflection at the Fermi energy, as one expects. Transparent objects are frictionless. 

To derive a q-swimming equation and fix the coefficients of Eq. ((21) we make use of the 
elementary observation that swimming is dual to pumping. A turning screw can be used to 
either pump or swim. The difference lies in the setup: In a pump the external forces and 
torques adjust to satisfy the constraints that the position and orientation of the pump are 
fixed while in a swimmer the position and orientation of the swimmer adjust to satisfy the 
constraint that there are no external forces and torques. 

Assume that no external forces are applied on a swimmer which can control its scattering 
matrix. The rate of momentum transfer is still given by Eq. Now, however, the energy 
shift has two terms: S = Sx + ^o- The first Sx (given in Eq. (jH))) arises from the swimmer's 
change of location, while the second Sq = iHSoSq comes from the swimming stroke. 

The total force acting on an adiabatic swimmer must vanish (to first order). This means 

n 

that V vanishes and the equation of motion for q-swimmers|15l|. is: 

= dEp\E) Trp (SoP) , (8) 

where t] is the friction coefficient, given in Eq. (|7j) and P is given in Eq. In the two 
channels case the trace on the right hand side is given by 

TrE{£oP) = hp{E) |r|2lm(rftlog(r/r')) . (9) 

We shall now apply the q-swimming equation, Eq. (jH)), to swimming in a Fermi gas at 
zero temperature in one dimension. This case is both simple and remarkable for, as we shall 
see, q-swimming turns out to be topological. A small deformation of the swimming stroke 
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does not affect the swimming distance which is an integer multiple of the Fermi wave length. 
This follows immediately from Eqs. ()8I9|) which combine to give: 

dX = ^lm id log(r/r')) =^da (10) 

where Ai? is the Fermi wave length. The fact that righthand side is an exact differential of 
the parameter a has two consequences: First, to swim one must encircle the point where 
the scatterer is transparent, r = r' = 0, and second, the distance covered in a stroke is 
quantized as a multiple of Xf/'^- The result is general and does not depend on the specifics 
of the swimmers. 

A swimmer will normally not have direct control on parameters a. Rather, it will control 
some physical parameters that will determine the scattering matrix (see examples below). 
Consider a swimmer with two independent controls. A stroke is a closed path in the plane 
of controls Since the reflection r is a complex valued function of the controls, one 

expects that by adjusting two controls, one can find points where r = . With each such 
point of transparency one can associate an (integer) index that counts how many times r 
rotates around the origin in one cycle around the point. We call the index the vorticity. The 
swimming distance in a closed stroke is proportional to the vorticity enclosed by the path. 

Let us now consider two examples of topological swimmers. 

Pushmepullyou: Consider 'a molecule' made of two scatterers. The scattering matrices 
associated with each scatterer have aj = 7j = for both scatterers and ri = cos^i, r2 = 
cos^^2- The two scatterers are separated by distance £, see Fig. ([21). The swimmer can control 

and the ratio Txjr^- The total scattering matrix of the swimmer can then be computed 
by considering the multiple scattering processes between the two scatterers. A computation 
yields for the zeros of total reflection of r the solutions of: 

e^^^Vi - r2 = 0. (11) 

It follows that the vortices occur when ri = ±r2 and the distance 11 is an integer (half and 
integer) multiple of wavelengths. Since the zeros are simple the vorticities are ±1. 

Three linked spheres: Consider a swimmer (proposed by Najafi and Golestanian in 
in a different context), made of three identical scatterers on a line, separated by distances 
£i and £2, see Fig. |21 Take = 7^- = and = cos6', independent of j = 1,2,3. The 
swimmer, a vibrating "trimer", now controls only the two distances ^1,^2- A computation 
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ri/r2 

FIG. 2: The vortex structure, atT = 0, of PushmepuUyou (left) and the three linked spheres (right). 
The swimming, in both cases, is topological: The distance covered in one stroke is proportional to 
the sum of the vorticities encircled by the path (the ellipses in the figures). The red dots have 
vorticity +1 and the blue dots have vorticity —1. 

shows that the total reflection r vanishes provided 




FIG. 3: The curvature for the three linked spheres, shown in a 3D plot as a function of control 
parameters at low temperature, T = To/4, (above) and high temperature, T = 4To, (below). The 
height of the peaks at the lower temperature is about a million times the height of the peaks at the 
higher temperature. The curvature vanishes as the temperature gets higher. The distance covered 
in one stroke is the integral of the curvature over the surface in the space of controls enclosed by 
the stroke. 



In geometrical language, Eqs. ()8I10|) define a connection dX in the space of controls. The 
distance covered by the swimmer in a one stroke C is then given by AX = j^dX. When 
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FIG. 4: The curvature at finite temperatures for PushmepuUyou, shown in a 3D plot as a function 
of control parameters. The upper picture corresponds to the case of low temperature. The case 
of high temperature is shown in the lower picture. In contrasts with the three linked spheres, for 
PushmepuUyou the curvature does not vanish at high temperatures. 

there are only two control parameters {x,y) (as in the examples above), an application 
of the Stokes formula gives for the displacement AX = jj f {x , y) dx dy , where f{x,y) is 
the (scalar) curvature, and the domain of integration has the boundary C. A plot of the 
curvature / is often an efficient way to describe a swimmer, see e.g.. Figs. 13 141 

It is instructive to see how topological Fermi swimmers are affected by temperature. At 
temperature T a region proportional to T near the Fermi energy will contribute to the inte- 
gral in Eq. (jH)). This makes the curvature a smooth function (rather than a collection of delta 
functions) on the space of controls. The total curvature enclosed in a path will now depend 
smoothly on the path. The temperature scale is determined by Tq = h"^ /AmkB^^F^ where ks 
is Boltzmann constant and m the mass of the scattered particle. At this temperature the 
support of / near a vortex becomes comparable with the distance between vortices. From 
the definition of Tq it follows that the larger the swimmer (the larger i) the more sensitive 
it is to temperature. 

PushmepuUyou and the three linked spheres have an essentially different behavior at 
high temperatures. Since the neighboring vortexes of three linked spheres (at T = 0) are 
of opposite sign, the smearing at high temperatures leads to vanishing curvature at high 
temperatures, see Fig. IHl Thus the three linked spheres do not swim effectively at high 
temperatures. On the other hand, the smearing of the the vortices of PushmepuUyou (at 
T = 0) does not lead to mutual cancellation, see Fig. |3J That means PushmepuUyou can 
swim effectively also at high temperatures. 

In the course of its motion a swimmer will, in general, transfer energy to the medium 
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and dissipate energy. As we shall now show, swimming without dissipation is possible in 
one dimensional Fermi gas at zero temperature. 

We define the dissipation, D, as the difference between the outgoing and incoming energy 
current. From 

D = ^TrE,{SS*)>0. (13) 

One of the consequences of this relation is that a non-dissipating swimmer is indistinguish- 
able from a static scatterer (since S = 0). Clearly, there is no dissipation in the former, and 
there should therefore be no dissipation in the latter. 

To see why it is possible to swim without dissipation, write Eq. (|13|). (in the two channel 
case), in the form: 

D = ^{{a-2kpXf cos^9 + e^ + (14) 

The first term vanishes by Eq. (fTUj) . Thus if 6* = 7 = there is no dissipation. We call a 
swimming without dissipation "super-swimming". 

A super-swimmer needs a larger space of controls to ensure that 6' = 7 = 0. (For instance, 
to make a super-swimmer out of Pushmepullyou one needs to control also the overall phase 
e*"^ of the two individual scatterers.) Super-swimmers are, in general, not transparent and 
one still needs to invest power to drag them. Only when the super-swimmer swims on its 
own, there is no transfer of energy to the ambient medium. 

We note that no-dissipation, S = 0, implies Eq. (fTUI) and hence implies quantization. 
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